® 



J 



Eur p§isches7gt ntamt 
European Patent Office 
Office europ^n des brevets 



® Publication number: 



0 222 450 

A2 



® 



EUROPEAN PATENT APPLICATION 



© Application number: 86201940«3 


® lnt.Ci/:G01S 15/02 , G01S 7/62 


@ Date Of filing: 07.11,86 




(S) Priority: 07.11.85 US 795801 


® Applicant: North American Pfililps 


@ Date Of publication of application: 


Corporation 


100 East 42nd Street 9th Roor 


20.05.87 Bulletin 87/21 


New York. N.Y. 10017(US) 


@ Designated Contracting States: 


® Inventor: Barnes, Casper William 


DEFR QB 


INT. OCTROOIBUREAU B.V. Prof. Hoistlaan 6 




NL-5658 AA Elndhoven(NL) 




® Representative: Veenstra, Qustaaf et al 




INTERNATIONAAL OCTROOIBUREAU &V. 




Prof. Hoistlaan 6 




NL-5856 AA Etndhoven(NL) 



0 Method and apparatus for estimating the ultrasound attenuation • VS - frequency slope of a 
propagation medium from the complex envelope of a signal. 

® An ultrasound attenuation scanner measures the local attenuation of tissues from pulse echo A-line data. An 
adaptive filter (60) with a transfer function which is the inverse of the tissue (25) undergoing examination 
operates on the complex envelope (*o{n)) of the RF signal. The filter (60) together with an FM detector (65) and 
an accumulator (70) functions as a negative feedback loop which continually maintains the filter characteristic as 
the inverse of the frequency dependent characteristic of the tissue (25). 
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Method and apparatus f r estimating the ultrasound attenuation -vs -frequency slope of a propagation 

nnedium from the complex envelop f a signal. 

The inv ntion relates to a method for measuring the local slope of the ultrasound attenuation versus 
frequency curve of a propagation medium comprising the steps of: 
transmitting a pulse of radiofrequency ultrasound energy into the medium; 

detecting echoes of said pulse which are reflected from points within said medium and producing an 
5 electrical signal from said echoes. 

The invention also relates to an apparatus for measuring the local slope of the ultrasound attenuation 
versus frequency curve at points In an object comprising: 

means for transmitting pulses of radio frequency ultrasound energy into an object, for detecting echoes of 
said pulses which are reflected from points within the object and for producing an RF electrical signal which 
70 represents said echoes. The method and apparatus are particulariy useful for the implementation of 
ultrasound attenuation scanners for the diagnostic imaging of human and animal tissues and body . 
structures. 

Most echo ultrasound Imaging apparatus of the prior art generates images of structures within a body 
wherein the brightness of Individual pixels con-esponds to the magnitude of local discontinuities in the 
16 acoustic impedance at corresponding points in the body. 

A new class of ultrasound imaging equipment, which porudces images of structures within a body 
wherein the brightness of pixels corresponds to a local value of ultrasound attenuation at a corresponding 
point in the body have recently been described. U.S. Patent 4,389,893 is incorporated herein, by reference, 
as background material. That patent describes apparatus which utilizes differential measurements of echo 
20 intensities along A-lines to estimate local tissue attenuation. A raster scan is used to generate images from 
such A-line measurements. U.S. Patent 4.515,163 also describes a method for determining local tissue 
attenuation from the positions of zero crossings in A-iine signals. 

It is possible to directly determine the slope of the attentuation vs. frequency characteristic of a 
propagation m^jium from samples of a signal which has propagated through that medium. However, in the 
25 context of an ultrasound attenuation scanner, that method requires that a radio frequency ultrasound, signal, 
which has propagated through the medium, be sampled at a rate which is at least twice as high as its 
highest frequency component in order to avoid aliasing error. This requires high speed sampling circuits 
which can significantly increase the cost and complexity of a scanner. 

It is an object of the invention to provide a method of the kind set forth that can be implemented with 
30 simplified signal processing circuitry. It is a further object of the invention to provide an apparatus of the 
kind set forth that is suitable for [>erfonming this method. 

To achieve these objects the method according to the invention is characterized in that it further 
comprises the steps of: 

extracting the complex envelope of said electrical signal; 
36 filtering the complex envelope using a complex filter with a controllable transfer function which can 
approximate the inverse of the transfer function of the propagation medium; 

measuring the average frequency of the filtered signal, the value of said average frequency being a 
representation of the value of the said attenuation characteristic of the medium; 

integrating the value of the measured average frequency and applying said integrated value, as negative 
40 feedback to control the transfer function of the complex filter, whereby said transfer function is maintained 
as the inverse of the transfer function of the medium. 

The apparatus according to the Invention is characterized in that it further comprises: 
complex demodulator means which extract the complex envelope of the RF signal; 

complex compensating filter means, having a controllable transfer function which can approximate the 
45 inverse of the transfer function of the object, connected to filter the complex envelope; 

FM detector means connected to detect frequency deviations of the filtered complex envelope signals, said 
frequency deviations being a measure of tiie slope of said attenuation versus frequency curve; 
accumulator means for integrating said frequency deviations and for applying the integrated deviations as 
negative feedback to control tiie transfer function of the complex filter and maintain it as the inverse of the 
60 transfer function of tfie object. 

The method and apparatus of the invention operate directly on tine complex envelope of a radio 
frequency pulse-echo signal. The complex envelope can be sampled at a much lower rate than tfie RF 
signal wittiout suffering from aliasing error. The signal processing circuitry of tiie scanner is thus 
significantiy Amplified. 
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In the absence of frequency dependent attenuation, the average frequency of the complex envelope of 
an ultrasound pulse which propagates through a medium is zero. The presence of frequency-<j0pend nt 
att nuation causes the averag frequency of the complex envelope to shift away from zero. This variation Is 
measured using an FM detector. The output of the detector is used as negative feedback to control th 
transfer function a compensating filter so that the filter drives th av rage frequency of the shifted complex 
envelope back to zero. The feedl>ack signal thus sets parameters of the compensating filter. Those 
parameters are a measure of the attenuation parameters of the propagation medium. In a pulse echo 
system, the tissue characteristics change with depth and the para meters of compensating filter track them. 
The invention thus provides an system that continuously adapts itself to the properties of the medium 
undergoing measurement. 

The invention will now be explained in detail with reference to the drawings, in which: 
Rgures la. lb and 1c illustrate the power spectra of signals; 
Rgures 2a and 2b illustrate the principal of the compensating filter: 
Figure 3 illustrates digital filter design parameters: and 

Rgure 4 Is a block diagram of an ultrasound attenuation scanner which incorporates the present invention. 

A real, stationary random signal. v,{t), propagates through a medium in which the attenuation is a 
function of frequency. The resulting output signal is Vo(t). 

If the medium is linear and time-invariant, then its effect on the signal can be described by a 
convolution in the form: Vo(t) = fdJ g(t-T) Vi(T), (1) 

where g(t) is the impulse response of the medium. The transfer function of the medium is: 
G(f) = /dt e-'irift g(t). (2) 

The invention Infers properties of the medium attenuation from signal processing operations on Vi(t) and 

Vo(t). 

The complex analytic representation of the real signal. v(t), is given by: 
T(t) = v(t) + i ^ <v(t)>. (3) 

where ^ denotes the Hilbert transform. In polar form this equation can be expressed as 
y(i) = a(t)e'^(^)- (4) 

At a reference frequency, fo. 
♦(t) = fl(t)-2irfot. (5) 

where 

y(X) « a<t)e'*<Virifot. (6) 

The complex envelope of the signal is defined as ^(t) A a(t)e*^W- (7) 
the reference frequency is chosen so that 
«o = E { I (t)} (8) 
then, 

E{0(t))»O. ^ ^ 

In these equations E { q (t)} and E { are the expected values of the time derivatives of ${t) and 
♦(t) respectively. From these equations it is clear that the average slope of the signal ^(t) can be made zero 
by a proper choice of the reference frequency fo. 

The relationship between these signals is illustrated in the drawings where Rgure la represents th 
power spectrum of a real signal v(t), Rgure lb represents the power spectrum of the complex analytic 
representation ^(t) of the signal of Rgure 1a. and Rgure 1c represents the power spectrum of the compi x 
envelope (t) of the signal of Rgure la- 



Effect of the Medium on the Complex Envelope. 

Let -rrft) « Vi(t) + i ^<V| (t)> = UW-wifoX (9) 
Then v,(t)= Re {7i(t)} . (10) and 
V o(t) = /dT g(t-T) v,(T) 
= /<^T g(t-T) Re {*,(T)e'ii ifoT} (11) 
= Re {/dT g(t-T)e"»irifo(t-T) ^,{T)e',rifot} 
- Re {*o(t)e'irif ot} (12) 
Where 

*o(T) = /dT g(t-T)e-'irifo(t-T) ^,(T) (13) 
or ^o(T) = /dT k(t-T) i^,(T) (14) 
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Where \ 
m = 9(t) e-'irifot (15) 
K(f) = Q(f +fo) (16) 

5 

Compensating fitter 

The effect of linear frequency dependence of attenuation on the complex envelope of a signal can be 
cancelled by a complex filter which filters the complex envelope. 
10 If the transfer function of the medium is 
Q(f) A e'*^. f O, (17) 

then the transfer function of the compensating filter for the complex envelope has the form 
H(f) = t ek(f+fo).f>-fo(18) 

(see Rgure 2a, in which 1 is the medium, 3 is a complex demodulator and 5 is a compensating filter) or H- 
(f) = 

1 



20 



25 



30 



35 



(See Rgure 2b in which 7 is the medium and 9 is a compensator). 

A digital implementation of a compensating filter for the medium can be determined as follows: 
H(v) = BK(^+Vo)for-i<v<i(19) 
Where Vo -^-^ (^). 
fs= sample rafe (21), and 

= kf». (22) 

The transfer function of the filter is illustrated in Rgure 3. 
The complex impulse response of the filter is: 



li(ii) =/dv H(v) e TTinv (23) 



1/2 

^ w r V2 ^ 

11(11) = ©^o I dv o(k+2 H in)v (25) 
J -1/2 

« » ^ 

5 + 2 iTin (26) 

ix(n) = 2(-l)'-^e^o Six^h (k^2) ^^^j 
k + 2 ITin 

ix(n) = 2(-l)" Sinn (k/2) (^ _ 2-nln) (28) 

k2 ^ 4Tr2^2 
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Rgure 4 illustrates an implementation of an attenuation scanner. Electrical pulses from an ultrasound 
pulse transmitter 10 are directed through a TR switch 15 to an ultrasound transducer 20. The transducer 
produces puis s of ultrasound nergy which propagate Into an object 25 where they produce echoes which 
are reflected back to the transducer. The echoes produce electrical signals in the transducer 20 which are 
again directed through the TR switch 15 to a receiver 30 where they are amplified In a conventional 
manner. The echo signals may be processed in a conventional AM detector -TGC circuit 35 to produce an 
image which is presented on a CRT display 40 in a conventional manner. The transmitter, detector, and 
display operate under the influence of a scanning and control circuit 45 which also elecfrically or 
mechanically scans the beam of ultrasound energy produced by the transducer over a sector of the object 
25. 

In accordance with the present invention the output signal Vo(t) of the receiver 30 is also fed to a 
complex demodulator circuit 50 which, at its output, produces a signal which corresponds to the complex 
envelope of the received signal. Digital implementations of a complex demodulator are described, for 
example, in United States Patent 4,481,646 and in "Theory and Implementation of the discrete Hilbert 
Transform". Gold et al in Digital Signal Processing, L Rabiner. Ed; IEEE Press 1972. The output of the 
complex demodulator 50 is sampled, in a sampling circuit 55 to produce a signal ^o(n). The sampling circuit 
55 need only operate at a rate equal to twice the highest frequency component In the complex envelope. 
The output of the sampler 55 is fed to the input of a digital compensating filter 60, which is implemented in 
the manner described in the previous paragraphs. The output of the compensating fiter. iAi(n), represents 
the compensated complex envelope of the echo signal. If the digital compensating filter parameters are 
correctly chosen the compensated complex envelope will con^espond to the complex envelope of the 
transmitted signal v i (t) and will have a zero average frequency. The compensated complex envelope signal 
is fed to a digital FM detector 85. The output of the digital FM detector. ^ ,(n), corresponds directly to the 
local value of the slope of the frequency dependent attenuation of the propagation medium 25. it may be 
applied directly to the display 40, for example to modulate the color of the display produced by the 
conventional AM detector 35. 

The output of the digital FM detector 65 is also fed to the input of an accumulator circuit 70 having a 
transfer function 



1 



1- - jz _ 

The output of the accumulator ^c(n) is then applied, as a negative feedback parameter ^ « -i9*c(n) to 
control the characteristic K ■ of the digital compensating filter 60. If k < k t^en ♦c(n) will increase, 
causing ^ to decrease (overcompensated case). If > ^ then (n) will decrease, causing ^ to 
increase (undercompensated case). The filter 60. FM detector 65 and accumulator 70 thus function as a 
negative feedback loop which continually forces the average frequency of the output of the digital 
compensating filter to zero to maintain the filter characteristic as the inverse of the frequency dependent 
characteristic of the propagating medium 25. 



Claims 

1 . A method for measuring the local slope of the ultrasound attenuation versus frequency curve of a 
propagation medium comprising the steps of: 

transmitting a pulse of radio frequency ultrasound energy into the medium; 

detecting echoes of said pulse which are reflected from points within said medium and producing an 
electrical signal from said echoes; 

characterized in that the method further comprises the steps of: 
extracting the complex envelope of said electrical signal; 

filtering the complex envelope using a complex filter with a controllable transfer function which can 
approximate the inverse of the transfer function of the propagation medium; 

measuring the average frequency of the filtered signal, the value of said average frequency being a 
representation of the valu of the said attenuati n characteristic of th medium; 

integrating the value of the measured average frequency and applying said integrated value, as negative 
feedback to control the transfer function of the complex filter, whereby said transfer function is maintained 
as the inverse of the transfer function of the medium. 
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2. A method as claimed in Claim 1, characterized in that It further comprises the step of displaying the ^ 
value of the frequency deviation at pixels of an image whose coordinates are mapped to correspond to the 
locations of the reflecting points in the propagation medium. 

3. An apparatus for measuring the local slope of the ultrasound attenuation versus frequency curve at 
5 points in an object comprising: 

means for transmitting pulses of radio frequency ultrasound energy into an object, for detecting echoes of 
said pulses which are reflected from points within the object and for producing an RF electrical signal which 
represents said echoes, 

characterized In that the apparatus further comprises: 
w complex demodulator means which extract the complex envelope of the RF signal; 

complex compensating filter means, having a controllable transfer function which can approximate the 

inverse of the transfer function of the object, connected to filter the complex envelope; 

FM detector means connected to detect frequency deviations of the filtered complex envelope signals, said 

frequency deviations being a measure of the slope of said attenuation versus frequency curve; 
T5 accumulator means for integrating said frequency deviations and for applying the integrated deviations as 

negative feedback to control the transfer function of the complex filter and maintain it as the inverse of the 

transfer function of the object. 

4. An apparatus as claimed in claim 3. characterized in that it further comprises display means 
connected to present pixels in a displayed image which represent the value of the frequency deviations of 

20 the complex envelope and means which control the location of pixels in the displayed image to correspond 
with a mapping of the location of the associated reflecting points in the object. 

5. An apparatus as claimed in claim 3 or 4, characterized In that the transfer function, H(v). of the 
compensating filter means is: 



H(v) = e^(^-^^o) - i <^< i , 

where v© » jj-* , 
f 8 B a sample rite. 

fo = the nominal frequency of the RF signal, 
K = kfs* and 

k is an attenuation constant of the medium. 
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® Method and apparatus for estimating the ultrasound attenuation - VS - frequency slope of a 
propagation medium from the complex envelope of a signal. 



® An ultrasound attenuation scanner measures the 
local attenuation of tissues from pulse echo A-Iine 
data. An adaptive filter (60) with a transfer function 
which is the inverse of the tissue (25) undergoing 
examination operates on the complex envelope (yl^o- 
(n)) of the RF signal. The filter (60) together with an 
FM detector (65) and an accumulator (70) functions 
as a negative feedback loop which continually main- 
tains the filter characteristic as the inverse of the 
frequency dependent characteristic of the tissue 
(25). 



CO 
< 

CM 
CM 
CM 




FIG. 4 



LU 



Xerox Copy Centre 



0Vti European Patent 



Office 



EUROPEAN SEARCH REPORT 



Application number 
EP 86 20 1940 



DOCUMENTS CONSIDERED TO BE RELEVANT 



Category 



Citation ol document with indication, where appropriate, 
of relevant c 



Relevant 
to claim 



CLASSIFICATION OF THE 
APPUCATION (Int. CI.-*) 



JS-A-4 197 750 (HASSLER) 

Abstract; column 1, line 36 - 
:olumn 3, line 45; column 4, line 
LO - column 7, line 57 * 



JS-A-3 805 596 (KLAHR) 
^ Abstract; column 19, line 20 
column 20, line 63; column 24, 
Line 13 - column 31, line 33 * 



JS-A-4 043 181 (NIGAM) 

^ Abstract; column 6, line 

column 8, line 51 * 



23 



JS-A-4 543 82 6 (FERRARI) 

Abstract; column 3, line 39 - 
::olumn 5, line 61 * 



5'R-A-2 556 845 (CGR ULTRASONIC) 
* Abstract; page 5, line 3 - page 
7, line 3; page 7, line 18 - page 
L8, line 17 * 



SP-A-0 080 544 (IBM FRANCE) 
^ Abstract; page 5, line 21 
oage 7, line 21 * 6c US-A-4 481 
S46 (Cat. D) 

-/- 



The present search report has been drawn up for all claims 



1,3 



1,3 



1,3 



1,3 



1,3 



1,3 



G 01 S 15/02 
G 01 S 7/52 
G 01 S 7/62 



TECHNICAL FIELDS 
SEARCHED (Int. CI.4, 



G 01 S 

G 01 N 

G 01 H 

H 03 H 

H 04 L 



Place of search 
THE HAGUE 



□ate o( completion of the search 
22-07-1987 



Examiner 
OLDROYD D.L. 



CATEGORY OF CITED DOCUMENTS 

X : particularly retevant if tak n alone 

Y : particularly relevant if combined with another 

d cum nt of the same category 
A . technological background 
O : non- written disclosure 
P : mtermedi ate document 



T : theory or principle underlying the invention 
E : earlier patent document, but publish d on, or 

after the filing date 
D : document cited in the application 
L : document cited for other reasons 

& ; member of the same patent family, corresponding 
document 



European Patent 
Office 



EUROPEAN SEARCH REPORT 



DOCUMENTS CONSIDERED TO BE RELEVANT 



Applicatton number 
EP 86 20 1940 
Page 2 



Category 



Citation of document with indication, where appropriate, 
of relevant passages 



Relevant 
to claim 



CLASSIFICATION OF THE 
APPLICATION (im. CI.4) 



FUJITSU SCIENTIFIC AND TECHNICAL 
JOURNAL, vol. 21, no. 2, June 
1985, pages 165-181, Kawasaki, 
JP; H. MIWA et al . : "Attenuation 
extraction from 

quadrature-detected phase of an 
ultrasonic echo in bio-tissue 
characterization" 

* Page 165, abstract; pages 
167-172, sections 2-6 * 



IEEE TRANSACTIONS ON 
ULTRASONICS, FERROELECTRICS , AND 
FREQUENCY CONTROL, vol. UFFC-33, 
no. 3, May 1986, pages 280-286, 
IEEE, New York, US; G. 
SALOMONSSON et al . : "On the 
separation of attenuation and 
texture of tissue using a 
parametric time-varying network" 
* Page 280, abstract; pages 
283-285, section IV * 



US-A-4 202 215 (MEYER) 
* AlDstract; column 3, line 24 - 
column 4, line 53; column 5, line 
21 - column 16, line 21 * 



IEEE TRANSACTIONS ON SON ICS AND 
ULTRASONICS, vol. SU-21, no. 2, 
April 1974, pages 91-97, New 
York, US; D.H. McSHERRY: 
"Computer processing of 
diagnostic ultrasound data" 



1,3 



1,3 



TECHNICAL FIELDS 
SEARCHED (Int. 01.^, 



1-4 



The present search report has been drawn up for all claims 



Place of search 
THE HAGUE 



Date of completion of the search 
22-07-1987 



Examiner 
OLDROYD D.L- 



CATEGORY OF CITED DOCUMENTS 

X : particularly relevant if taken alone 

Y : particularly relevant if combined with another 

d cument of the same category 
A technological background 
O : non-written disclosure 
P intermediate document 



T : theory or principle underlying the invention 
E : earlier patent d cument. but published on. or 

after the filing date 
D : document cited in the applicatton 
L : document cited for other reasons 

& : member of the same patent family, corresp nding 
document 



European Patent 
Office 



EUROPEAN SEARCH REPORT 



DOCUMENTS CONSIDERED TO BE RELEVANT 



Application numDer 

EP 86 20 1940 
Page 3 



Category 



Cjtatton of document with indication, where appropriate, 
of relevant passages 



Relevant 
to Claim 



CLASSIFICATION OF THE 
APPLICATION {Int. a.4) 



SLECTRONICS LETTERS, vol. 8, no, 
L5, 27tli July 1972, pages 
380-382, Hitchin, GB; G.D. CAIN: 
'*Hilbert-t.ransf orm description 
linear filtering" 



JLTRASONIC IMAGING, vol. 4, 
1982, pages 355-377, Academic 
Press Inc., New York, US; K.V. 
^URUMURTHY et al. : "A dispersive 
itiodel for the propagation of 
iltraspund in soft tissue" 
Whole article * 



TECHNICAL RELDS 
SEARCHED Ont. O.^i 



The present search report has been drawn up for all claims 



Place of search 
THE HAGUE 



Date of completion of the search 

22-07-1987 



Examiner 

OLDROYD D.L. 



CATEGORY OF CITEO DOCUMENTS 

X : particularly relevant if taken alone 

V - particularly relevant if combined with another 

document of the same category 
A : technol gica) background 
O : n n-written disclosure 
P intermediate document 



T : theory or principle underlying the invention 
E : earlier patent document, but published on. or 

after the filing date 
D : document cited in the application 
L : document cited for other reasons 

& : member of the same patent family, corresponding 
document 



